Insulating behavior at the neutrality point in dual-gated, single-layer graphene 
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The fate of the low-temperature conductance at the charge- neutrality (Dirac) point in a single 
sheet of graphene is investigated down to 20 mK. The potential fluctuations in graphene are made 
artificially small by using a dual top- and back-gate geometry such that screening reduces the 
intrinsic disorder potential. As the temperature is lowered, the peak resistivity diverges with a 
power-law behavior and becomes as high as several Megaohms per square at the lowest temperature, 
in contrast with the commonly observed saturation of the conductivity. As a transverse magnetic 
field is applied, our device remains insulating and directly transitions to the v—0 quantum Hall state 
at a field of 100 mT. 



The ability to create electronic devices in graphene has 
made it possible to study 2D Dirac fermions in the solid 
state [TH3] . Transport measurements in a large magnetic 
field display quantum Hall plateaus with unconventional 
values of conductance, a signature of the Dirac equa- 
tion describing electrons in graphene pQ. Such behavior 
emerges as the cyclotron gap opens in a magnetic field: 
when this gap becomes larger than disorder-induced fluc- 
tuations in the surrounding potential, the effect of the 
linear Dirac band structure becomes evident. At zero 
magnetic field, the disorder landscape dominates [4) [30], 
blurring the interesting phenomena that might occur at 
the Dirac point. Recently, the influence of disorder has 
been reduced by either suspending a sheet [5 or placing 
it on atomically-flat boron nitride [6] , and many discover- 
ies in transport have been made due to the more readily- 
accessible Dirac point in these cleaner systems [7HTD]. 

The nature of the conductivity at the charge-neutrality 
point <jnp has been debated since the first graphene- 
based devices were fabricated. Theory for ballistic 
graphene predicts a value of 4e 2 /irh for cfnp [HI fl2] . 
However, early experiments on graphene measured a^p 
in the range of 2-12 e 2 /h [13] Q3]. It was soon re- 
alized that <jnp was sample-dependent and determined 
by the density of carriers in electron and hole puddles 
due to disorder from static charges on or near the sheet 
of graphene p~3j [15j [16]. In suspended graphene de- 
vices [15 the conductivity showed a more pronounced 
temperature dependence, but still saturated at low tem- 
perature and remained higher than 4e 2 /7rh. Recently it 
has been shown that the potential landscape in graphene 
can be made artificially clean by screening potential fluc- 
tuations with a second nearby, doped graphene sheet [17] . 
In that work, instead of saturating at values near e 2 /h, 
<jn p dropped with a power-law temperature dependence 
T a , where a=2 for the most insulating samples, down to 
T=4K. Further, the authors observed a strong magne- 
toresistance in the temperature regime above 10 K and 
attributed it to weak localization, inferring that ultra- 
clean graphene may be an Anderson insulator. However, 
a recent theory [18] postulates that this temperature de- 



pendence results from a "disorder by order" effect in 
graphene, where increasing order causes the sample to 
become more insulating at low density. Here a^p oc T a 
naturally emerges as the temperature dependence of con- 
duction through a landscape of electron and hole pud- 
dles. To distinguish between the two competing scenar- 
ios, lower-temperature transport, where the electronic co- 
herence important in localization is more prominent, is 
needed. 

In this Letter, we report on electronic transport in 
dual-gated, single-layer graphene devices, in which the 
back gate is used to control the global carrier density, 
and the top gate is primarily used to screen the Coulomb 
disorder potential, not to locally control the carrier den- 
sity. At the charge neutrality point (CNP), the de- 
pendence of the conductance qnp on temperature (T) 
and magnetic field (B) is investigated down to temper- 
atures of 20 mK. The temperature dependence is strong 
down to T=400mK and can be fitted to a power-law 
g NP oc T a with a^0.48±0.05. The low-field magne- 
toresistance shows an uncommon behavior around B=0 
that is inconsistent with weak localization. The quantum 
Hall regime is entered at just 0.1T, but instead of show- 
ing a conductance plateau at 2 e 2 /h around the CNP, the 
device is insulating, indicating that the spin and valley 
SU(4) symmetry is broken at very low fields. A peak 
in gNp(B) marks a transition between the zero field in- 
sulating state and the v=0 insulator at a field on order 
100 mT, about 20 times lower than for our non top-gated 
devices. From the low-temperature saturation of gNp(T) 
and the field at which the u=0 quantum Hall plateau 
emerges, we estimate that the potential fluctuations are 
less than 100 /ieV, an order of magnitude lower than re- 
ported for the suspended devices. 

We fabricated our devices using hexagonal-boron ni- 
tride (h-BN) as a substrate for graphene, with good 
electronic properties [6j [19] resulting from the extreme 
flatness and cleanliness of h-BN flakes. Single-layer 
graphene flakes were exfoliated on a stack of polyvinyl 
alcohol (PVA) and polymethylmethacrylate (PMMA). 
Separately, h-BN was exfoliated on a silicon wafer with 
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FIG. 1: (a) Schematic of the device in voltage biased mode. A 
bias vac is applied to the sample, the current %r> is then collected 
at the drain and measured with a lock-in amplifier. A voltage Vbg 
is applied to the degenerately doped piece of wafer to control the 
carrier density in the whole sheet. Vtg is applied to a suspended 
metallic gate, 70nm away from the graphene sheet, which varies 
the carrier density underneath it. (b) Conductivity a as a function 
of both gate voltages, measured at a temperature T=4K. 



a 300nm-thick thermal oxide, then exposed to flowing 
Ar/02 at 500 °C for 8 hours, which removes organic con- 
tamination from tape residue [2TJ [22]. Each flake was 
characterized by atomic force microscopy and Raman 
spectroscopy prior to transfer [20 . The PVA layer was 
dissolved in deionized water, which lifted off the PMMA 
membrane with the attached graphene flake [6 j. The 
graphene flake was then aligned to the h-BN, pressed 
on top and baked at 110 °C to promote adhesion be- 
tween the two flakes. The PMMA membrane was dis- 
solved in acetone/IPA and the graphene/h-BN stack was 
again treated in Ar/0 2 at 500 °C for 4 hours [21 . The 
cleanliness of the whole stack was confirmed with Raman 
spectroscopy [20 . Successive electron-beam lithography 
steps were used to oxygen-etch the graphene flake into 
a Hall-bar and to make electrical contact to the device 
(1 nm Cr, 200 nm Au). Finally, a 1 /im wide top gate, also 



patterned by electron-beam lithography, was suspended 
70 nm above the flake [23H26] . For a last cleaning step, 
we could not anneal the device in oxygen, as we found it 
damages the metallic contacts [21]. Instead, the device 
was annealed at 325 °C in Ar/H2 for 4 hours. Additional 
details on the fabrication and characterization are avail- 
able in Ref. [20]. A schematic of the completed device 
geometry is shown in Fig. 1(a). Another device, called 
the control device, was fabricated with no top-gate on 
the same sheet of graphene but electrically isolated from 
the top-gated device. The control device showed none of 
the distinct transport properties of the top-gated device 
at the CNP, helping discern which transport features are 
solely due to the presence of the suspended gate. 

Graphene devices were measured in two differ- 
ent cryostats: a variable-temperature insert enabling 
temperature-dependent transport measurements from 
300 K down to 1.7 K, and a dilution fridge where samples 
were measured at lower temperatures, down to 20 mK. 
The conductance g is determined in a standard voltage- 
biased lock-in measurement with an excitation voltage of 
4 /xV at 92.3Hz. The resistance r is defined as 1/g. DC 
voltages are applied to the top-gate (Vtg) and back-gate 
{Vbg). 

g is shown on Figure 1(b) as a function of Vtg and 
Vbg at T=4K. The carrier density can be controlled 
independently and with opposite polarities underneath 
and outside the top-gated region. As in previous work 
on dual-gated graphene [3j [271 EH] > 9 exhibits local min- 
ima along two intersecting lines corresponding to each 
region being tuned through the CNP. However, unlike 
in typical dual-gated graphene devices, gwp is near zero 
along these lines, g was also measured in a 4-probe ge- 
ometry, from which we extracted the resistance of the 
cryostat's wiring, 2.5 kfi. Underneath the top-gate, 
CtgVtg + CbgVbg — at the CNP, yielding a top- 
gate-to-back-gate capacitance ratio of 1.3 from the slope 
of the diagonal line in the (Vbg, Vtg) plane. Cbg is 
5.94(±0.5) x 10 10 cm _2 V _1 , as extracted from the period- 
icity of Shubnikov-de-Haas oscillations. Using a parallel- 
plate capacitance model, we estimate that the top gate 
is 70 nm away from the flake, which was confirmed by 
atomic force microscopy. The device exhibits little intrin- 
sic doping, with a CNP voltage of -2.2 V on the back-gate, 
and a mobility of 60,000 cm 2 /Vs, as extracted from a lin- 
ear fit to g(VsG) at the CNP. We note that the mobility 
of top-gated and non top-gated regions was comparable, 
which shows that the suspended gate does not deteriorate 
the electronic properties of our device. 

Unlike typical graphene samples, r as a function of 
Vbg has a strong temperature dependence in our device 
(all measurements from here on are taken with Vtg=0). 
The resistance at the CNP tnp dramatically increases 
at low temperature, from 13 k^ at T=300K to 300 k^ at 
T=4K [inset, Fig. 2(a)]. By contrast, the peak in r of 
the control device is only 10 k^ at 4K [20], comparable 
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FIG. 2: (a) Resistance as a function of Vbg a t T=300K (dashed 
grey line) and T=2K (black line). Insert: Resistance at the 
charge neutrality point (Vbg — —2.2V) as a function of tempera- 
ture, (b) Conductance as a function of temperature on a log scale, 
at Vbg — — 2.2V. Open circles and filled squares correspond to 
data points taken in the dilution fridge and variable temperature 
cryostat respectively. Error bars correspond to one standard devi- 
ation: for T>5K, these are smaller than the dots. The grey line 
corresponds to a fit cr NP oc T a , with a w 0.48 for T<80K. At 
higher temperature cr^p rises with a faster exponent. 

to what is commonly seen in good-quality, single-layer 
graphene devices, r^p was measured at lower temper- 
ature in a separate cool-down using a dilution fridge: 
further lowering T to 400 mK increases r^p, at which 
point it measures 2.3 M^, then remains constant down 
to 20 mK within experimental error. The full-width at 
half-maximum in vnpCVbg) a t 20 mK corresponds to a 
residual impurity density on order 10 10 cm -2 . For carrier 
densities higher than ~5xl0 10 cm -2 , r decreases when 
T is lowered, similarly to what has been observed else- 
where on high-quality samples without a top gate [T5] 
and in our control device [20]. Fig. 2(b) shows ^atp(T), 
which follows a power law g^p oc T a as a function of 
the temperature, with a=0.48±0.05, as extracted from 
a linear fit to the curve. This insulating behavior is not 
due to the opening of a band gap, which would lead to 
an exponentially-activated conductivity. The tempera- 
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FIG. 3: Conductance at the charge neutrality point g^p as a 
function of the magnetic field at two different temperatures, T=2K 
(dashed grey line) and T=400mK (black line) 

ture dependence is also slower than the T 2 dependence 
measured in Ref. [17 (and expected from the Boltzmann 
equation with electron-electron scattering) . An evolution 
similar to ours was reported in suspended graphene [15] , 
although the overall conductance was several orders of 
magnitude higher. 

g N p(B) for T=2K and 400 mK in the range 
| B | ^250 mT exhibits 2 strong peaks located symmetri- 
cally at 5=±70mT (Fig. 3). The peak conductance 
at ±70 mT has less T dependence than gNp(B=0). A 
strong magneto-conductance was seen in Ref. [17] as well 
and was fit to the weak localization theory for graphene. 
However, our observed magneto-conductance is distinct 
from that seen previously in Ref. [17 or in typical two- 
dimensional conductors in two ways. First the magneto- 
conductance around B=0 has positive curvature, differ- 
ent from the cusp-like, negative curvature typical of WL. 
Second, after the peak at ±70 mT, the conductance drops 
abruptly towards zero, indicating that the device returns 
to an insulating state. It appears that the magneto- 
conductance peak is not due to Anderson localization, 
but rather marks a transition between the zero-field in- 
sulating state and the v = insulator. 

To investigate this transition, a low field fan diagram 
[g(^BG 9 B)] is measured at T=2 K and shown in Fig. 4(a). 
Away from the CNP, we observe plateaus for v = 2, 6, 10 
(dashed lines) that are well-developed on the hole side for 
B>0.bT. The cut g(VeG 7 B=l T) shows these plateaus in 
addition to the v=0 plateau around the CNP [Fig. 4(c)]. 
Quantum Hall plateaus are better resolved on the hole 
side than on the electron side, which we found to be com- 
mon for two-terminal quantum Hall conductance mea- 
surements on other graphene on boron nitride devices. 
Interestingly, the broken-symmetry v=0 state seems to 
persist all the way down to very low fields (solid blue 
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FIG. 4: (a) Conductance g as a function of the transverse magnetic 
field B and VbGj a ^ T=2K, measured in a two probes configuration. 
The transitions between plateaus are underlined for v =0, 2, 6, 10. 
(b) Cut of g(B) at the charge neutrality point Vbg — — 1-1V. (c) 
Cut of g(V B o) at B=1T. 



region that runs between Vbg ^0 an d -IV). 

Discussion. Typically the disorder in the sample may 
be extracted from the temperature below which g(T) sat- 
urates [5 . The temperature dependence in our device is 
very strong down to 400 mK, at which point it may 
remain flat down to 20 mK, the dilution fridge's base tem- 
perature, although the temperature dependence in this 
range is harder to track due to current noise in our setup. 
The temperature dependence is expected to be flat when 
fc#T is below the Fermi energy's fluctuations 5Ep, allow- 
ing us to coarsely estimate SEp and find an unexpectedly 
low value 5Ep ~ 40 \ieV . 

Another way to estimate 5Ep is to look at the mag- 
netic field at which the v=0 plateau develops. For our 
control device [20] the two-terminal conductance devel- 
ops plateaus at a moderate magnetic field (^1-2 T) with 
value 2e 2 /h around the CNP (corresponding to v=±2). 
In typical graphene-on-BN devices, we and others [33] 
find that the two-probe conductance at the CNP only 
starts decreasing from 2e 2 /h to zero for 5^2-3 T 20 , 
as the valley degeneracy of the n=0 Landau level is 
lifted [30]. The v — gap is usually fully developed 
at around 5T, at which point the current becomes un- 
measurable with our setup. This is qualitatively dif- 
ferent from what we observe with our top-gated de- 



vice, where the conductance around the charge neutral- 
ity point is never quantized and always very small com- 
pared to 2e 2 /h [Fig. 4(b)]. The conductance is heav- 
ily suppressed on the high field side of the conductance 
peak discussed earlier, corresponding to an opening of 
the v = gap at ~100mT. In a magnetic field, the en- 

2 

ergy of Coulomb interactions is on the order Ec = J^-, 

where Ib is the magnetic length Ib = \J^>- While to first 
order these interactions preserve valley symmetry, it has 
been shown that higher-order terms break this symmetry 
and are on order SEc = ^-^c, where a is the spacing 
between neighboring carbon atoms [31-33 . A naive es- 
timate of this contribution gives 5Ec^ meV/T. Inter- 
estingly, the field dependence of the v=0 gap has been 
studied in Ref [33], where it was found that the effec- 
tive g factor g/\ = dA Q /dB for the v = state had 
this same order of magnitude. We can use this result to 
estimate the Landau level broadening due to the Fermi 
level fluctuations from the onset of the v = gap at 
100 mT. We find SEp^lOO /iieV , in good agreement with 
our previous estimate based on the low-temperature sat- 
uration of gNp(T). This confirms that the Fermi energy 
fluctuations are at least one order of magnitude lower 
in our double-gated device than they are in suspended 
graphene [15] [34]. 

Using Ep = UvF^fwri, our inferred fluctuations in Ep 
correspond to density fluctuations on the order of 10 6 
cm -2 , two orders of magnitude lower than in most re- 
ported suspended devices [34 ! The density of impuri- 
ties is unlikely to be this small, but probably on order 
10 9 to 10 10 cm -2 [6j [19] and similar to our control de- 
vice. This is internally consistent: while Coulomb in- 
teractions are usually barely screened at the CNP [35] . 
the nearby grounded metallic plane provides additional 
screening and the Fermi energy fluctuations are therefore 
expected to be weakened. 

The improved screening is likely to be responsible for 
the diverging resistivity we observe. As charge puddles 
get shallower and further apart, the minimum conductiv- 
ity is expected to drop as it becomes harder for electrons 
to percolate [I6j[36]. We note that cjnp was calculated 
in the percolating regime [36 and predicted to depend 
on density fluctuations as cjnp ^ Sn 0A1 , a result that 
may be relevant here as the resistance at p-n boundaries 
between charge puddles grows to dominate sample resis- 
tivity [37]. However, the top gate is tens of nanometers 
away from the flake, significantly further than reported in 
Ref. [17] , and it is therefore surprising that the screening 
should have such dramatic consequences, as interactions 
should naively be screened only on distances larger than 
the spacing to the top gate. Further work is therefore 
needed to better understand the role of screening in such 
systems. 
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SUPPLEMENTARY INFORMATION 
Device fabrication 

Polyvinyl alcohol (2% in water) is spun at 6000 rpm 
on bare silicon and baked at 160 °C for 5 minutes, result- 
ing in a 40nm thick layer. Then, a layer of PMMA (5% 
in anisole) is spun at 2400 rpm and baked 5 minutes at 
160 °C. The total polymer thickness is on the order of 
450nm. Graphene is then exfoliated on this stack using 
Nitto tape, located using optical microscopy, and char- 
acterized with Raman spectroscopy. Boron nitride flakes 
are exfoliated on a silicon wafer piece with a 300nm thick 
thermal oxide, then baked in flowing Ar/02 at 500 °C to 
remove tape residue. The cleanliness and thickness of the 
flakes are characterized with atomic force microscopy and 
Raman spectroscopy. 

The PVA layer is dissolved in deionized water at 90 °C, 
which lifts-off the PMMA membrane with the graphene 
attached to it. The membrane is then adhered across a 
hole in a glass slide [SI] and baked at 110 °C. We then use 
a home-made probe-station to align the graphene flake 
on top of the boron nitride substrate. Once both flakes 
are in contact, the stack is baked at 120 °C to promote 
adhesion. The PMMA layer is dissolved in hot acetone, 
then rinsed in IPA, which leaves the graphene flake on 
top of the boron nitride flake. This stack is annealed 
in flowing Ar/02 at 500 °C for 4 hours, which removes 
process residue and leaves the graphene flake pristine, as 
checked by Raman spectroscopy. 

We use regular electron beam lithography combined 
with oxygen plasma etching to pattern a graphene Hall 
bar. In order to fabricate a suspended top gate above the 
device, the samples are spin-coated at 6krpm with a so- 
lution of polymethylmethacrylate (950k), 3% in anisole, 
then baked at 160 °C for 5 minutes. An additional layer 
of methyl-methacrylate (8.5% in ethyl lactate) is spin- 
coated at 6krpm and baked 5 minutes at 160 °C. The 
MMA layer is 50% more sensitive to electron irradiation 
than the PMMA layer and it is therefore possible to de- 
velop the top resist layer without exposing the bottom 
layer. The e-beam writing system we used is a JEOL 
6300, with an acceleration voltage of lOOkeV. The con- 
tacts and the feet of the suspended bridge are exposed 
with 1000/iC/cm 2 , which is enough to dissolve both re- 
sist layers upon development in MIBK/IPA (1:3) for 45s. 
The span of the suspended bridge is only exposed with 
a base dose of 650/iC/cm 2 , which only develops the top 
resist layer. After development, the device is cleaned for 
2 minutes with UV ozone, then metallized with lnm of 
chromium and 200nm of gold. 

Figure 5 is a scanning electron micrograph of a 
graphene on boron-nitride device with a suspended top- 
gate that has been fabricated using the same recipe. Sus- 
pended gates as long as 7 microns have been fabricated 




FIG. 5: Scanning electron micrograph of a graphene on boron 
nitride device with a suspended top gate. Scale bar=l ^m. 
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FIG. 6: a) Raman spectrum of the boron nitride flake used for 
the device described in the main paper, prior to transferring the 
graphene flake. The boron nitride Raman peak is labelled Gbn- 
b) Raman spectrum of the whole graphene on boron nitride device 
after oxygen annealing and before contacts were made to the flake. 
The peaks labelled G and 2D are attributed to graphene. 



using this recipe. These usually resist further heat treat- 
ment as well as cryogenic temperatures. We found that 
gate voltages as high as 40V can be applied to the top- 
gate without damaging it. 
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FIG. 7: Resistivity as a function of the back-gate voltage for the 
control device 

Raman spectra 

We check the cleanliness of every boron nitride flake 
we use with atomic force microscopy and Raman spec- 
troscopy. Figure 6. a shows the Raman spectrum of the 
boron nitride flake used in the device studied in the pa- 
per, prior to transferring the graphene flake. The absence 
of a broad background signal [S2] is a very clear indica- 
tion that the flake is free of any kind of organic contam- 
ination. The Raman spectrum of the transferred device 
before the last lithographic step -when the suspended 
gate is patterned, is shown on Figure 6.b. The ratio of 
the amplitudes of the 2D and G peak is I2D Ag = 5.5, and 
the full half width of the 2D peak is 19 cm -1 , indicating 
with no ambiguity that the flake studied here is single- 
layer graphene [S3] . The absence of a broad background 
signal attests for the cleanliness of the device. 
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FIG. 8: a) Two terminal inductance as a function of the transverse 
magnetic field and the back-gate voltage, b) Cuts of the conduc- 
tance as a function of the back gate voltage at B=2T (black line) 
and B=4T (red line). 



Resistance of a non top-gated sample 

The resistivity of a non top-gated part of the sample 
described in the main paper is shown on figure 7, mea- 
sured at T=4K. The peak resistivity is ^lOkSl and does 
not diverge as the temperature is lowered. 

Quantum Hall conductance of a non top-gated 
device 

The two-terminal quantum Hall conductance of a 
graphene on boron nitride device with no top-gate is 
shown on Figure 8. The conductance is measured at 
T=2K in a two-probe configuration. Most of the mea- 
surements in this paper use this geometry, as voltage- 
biased measurements of the conductance are more con- 
venient than current-biased ones when the device is in- 



sulating, around the charge neutrality point. 

As the magnetic field increases, the degeneracy of 
each Landau level splits and we observe new quantized 
plateaux. In particular, the v=0 phase slowly appears 
around B=3T: this is in sharp contrast with the abrupt 
transition seen in figure 4b of the main paper, which oc- 
curs at a much lower field. Figure 8b shows two sections 
of the Fan diagram at B=2T and B=4T. The conduc- 
tance at 2T shows the standard sequence of plateaux, as 
seen in other two terminal devices with no splitting of the 
Landau levels [S4]. We stress that the small dip in con- 
ductance at the neutrality point at 2T is only an artifact 
from the two terminal geometry, as shown in Ref. [S4], 
and is not related to the opening of the v=0 gap. At 
4T, the device starts being insulating around the charge 
neutrality point, as the v=0 gap opens up, and the other 
plateaus become further resolved. 
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